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Entanglement does not correspond to any observable and its evaluation always corresponds to an 
estimation procedure where the amount of entanglement is inferred from the measurements of one 
or more proper observables. Here we address optimal estimation of entanglement in the framework 
of local quantum estimation theory and derive the optimal observable in terms of the symmetric log- 
arithmic derivative. We evaluate the quantum Fisher information and, in turn, the ultimate bound 
to precision for several families of bipartite states, either for qubits or continuous variable systems, 
and for different measures of entanglement. We found that for discrete variables, entanglement may 
be efficiently estimated when it is large, whereas the estimation of weakly entangled states is an 
inherently inefficient procedure. For continuous variable Gaussian systems the effectiveness of en- 
tanglement estimation strongly depends on the chosen entanglement measure. Our analysis makes 
an important point of principle and may be relevant in the design of quantum information protocols 
based on the entanglement content of quantum states. 

PACS numbers: 03.67.Mn, 03.65.Ta 
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I. INTRODUCTION 

Entanglement is perhaps the most distinctive feature 
of quantum mechanics, and definitely the most relevant 
resource for quantum information processing [T'|. Indeed, 
quantification of entanglement and schemes for its mea- 
surement have been the subject of extensive efforts in 
the last decade 11 ii, S i, 0, S i, 0, [O, [H E Q 
m, [M EE m, HliJ Indeed, the entanglement con- 
tent of a quantum state is a crucial piece of information 
in the design of quantum information protocols, and a 
question naturally arises on whether quantum mechanics 
itself poses limits to the precision of its determination. As 
a matter of fact, any quantitative measure of entangle- 
ment corresponds to a nonlinear function of the density 
operator and thus cannot be associated to a quantum 
observable. As a consequence, any procedure aimed to 
evaluate the amount of entanglement of a quantum state 
is ultimately a parameter estimation problem, where the 
value of entanglement is indirectly inferred from the mea- 
surement of one or more proper observables. An opti- 
mization problem thus naturally arises, which may be 
properly addressed in the framework of quantum estima- 
tion theory (QET) |2l|, which provides analytical tools 
to find the optimal measurement according to some given 
criterion. 

Our aim is indeed to evaluate the ultimate bounds to 
precision posed by quantum mechanics, i.e the smallest 
value of the entanglement that can be discriminated, and 
to determine the optimal measurements achieving those 
bounds. Being entanglement an intrinsic property of 
quantum states we adopt local quantum estimation the- 
ory, where the optimal estimators arc those maximizing 
the Fisher information [22, 23, 24] and in turn minimizing 
the variance at fixed value of entanglement. Local QET 
provides any family of quantum states with a geomet- 



ric structure based on distinguishability and accounts for 
the optimal measurement that can be performed on the 
quantum system as well as the optimal data processing 
of the outcomes of the measurement. 

Local QET has been applied to the estimation of quan- 
tum optical phase [25j as well as to estimation problems 
involving non unitary processes in open quantum systems 
[26| . either in finite dimensional systems [13] or contin- 
uous variable ones [2^. This includes the optimal es- 
timation of the noiseparameter of depolarizing [29^ or 
amplitude-damping [2^ [s^] channels. Recently, the ge- 
ometric structure induced by the Fisher information it- 
self has been exploited to give a quantitative operational 
interpretation for multipartite entanglement [3l| and to 
assess quantum criticality as a resource for quantum es- 
timation [sl ]. 

In this paper we systematically apply local QET to 
the problem of efficiently estimate the amount of entan- 
glement of a quantum state. We consider several fam- 
ilies of bipartite states, either for qubits or Gaussian 
states, and evaluate the symmetric logarithmic derivative 
to estimate entanglement through different measures, e.g. 
negativity or linear entropy. Then we explicitly calculate 
the quantum Fisher information and derive the ultimate 
bounds to the precision of estimation. Overall, we found 
that, both for qubits and Gaussian states, entanglement 
may be efficiently estimated when it is large. On the 
other hand the estimation of a small amount of entan- 
glement for qubits is an inherently inefficient procedure, 
i. e the signal-to-noise ratio is vanishing for vanishing en- 
tanglement, whereas for Gaussian states it depends on 
the chosen measure of entanglement. We also found that 
the presence of other free parameters besides entangle- 
ment does not generally influence the estimation preci- 
sion, thus preventing the possibility of further optimizing 
the estimation procedure. 
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The paper is structured as follows: in the next Section 
we give some basic elements of quantum estimation the- 
ory and introduce the quantum signal-to-noise to assess 
the estimability of a parameter. In Section lTTTl we analyze 
the estimation of entanglement by means of negativity 
[33I I and linear entropy for the family of pure two-qubit 
states as well as for two families of entangled mixtures. In 
Sectio n- Hyl we address a family of PPT bound-entangled 
states [3JI for two-qutrit systems as an example of states 
with an inherently small amount of entanglement. In Sec- 
tion |V] we address entanglement estimation for Gaussian 
states, either pure states (twin-beams) or entangled mix- 
tures. Section IVll closes the paper with some concluding 
remarks. 



and, in turn, that Var(A) > [MH{X)]~^ represents the 
quantum version of the Cramer- Rao theorem, i.e. the 
ultimate bound to precision for any quantum measure- 
ment aimed to estimate the parameter A. The SLD itself 
provides an optimal measurement, that is, using a mea- 
surement described by the projectors over the eigenbasis 
of Lx we saturate inequality ([5|). 

Upon diagonalizing gx = J2kPk\^k){i'k\ and using 
Eqs. ([3|) and ([5]) we obtain 



r O {'4'n\dxgx\lpm/ I , w , 

La = 2 2^ — I — — IVn) Wr, 



Pn 



(6) 



Pn 



Pn + Pni 



II. QUANTUM ESTIMATION THEORY 

In an estimation problem one tries to infer the value 
the of a parameter A by measuring a different quantity 
X, which is somehow related to A. An estimator A = 
A(a;i, a:2, ..) for A is a real function of the outcomes of 
measurement. The Cramer- Rao theorem [s^ establishes 
a lower bound for the variance Var(A) of any unbiased 
estimator 

in terms of the number of measurements M and the so- 
called Fisher Information (PI) 



i^(A) = ^p(x|A) [dx\npix\\)r ^ 



[dxp{x\XW 
pix\X) 



(2) 



where p{x\X) denotes the conditional probability of ob- 
taining the value x when the parameter has the value 
A. 

In quantum mechanics, according to the Born rule we 
have p{x\X) = Ti[Exgx] where {E^} are the elements 
of a positive operator-valued measure (POVM) and gx 
is the density operator parametrized by the quantity we 
want to estimate. Introducing the Symmetric Logarith- 
mic Derivative (SLD) Lx as the operator satisfying the 
equation 



Lxgx + QxLx 



dgx 
dX 



(3) 



we have that dxp{x\X) — TrldxgxEx] — Re(TT[gxLxEx]), 
and the Pisher Information in Eq. Q may be rewritten 



as 



Re{TT[gxLxE,]f 

TiigxE^] 



(4) 



Starting from Eq. ^ one may prove that F{X) is upper 
bounded by the so-called Quantum Fisher Information 
(QPI) [2l,[2| 



Fix) < HiX) = TiigxLj 



(5) 



which, for a family of pure states gx 
to 



I'/'a)('0a| reduces 



Lx = 2dx\^x){^x\ 
H{X)^4[{dxtPx\dxi'x) 



(7) 



(dx^Pxlipx 



When more than a parameter are involved we have quan- 
tum states PA depending on a set of N parameters 
A — {Xj}, j = 1, . . . , N. In this case the geometry of 
the estimation problem is contained in the QPI matrix, 
whose elements are defined as H{A)ij — ^Ti[gA{Li, Lj}] 
where Li is the SLD that corresponds to the parameter 
Xi and {A, B} — AB + BA denotes anti-commutator. 
The explicit formula for the QI matrix reads as follows 



{diPn)(d]Pn) 



E 



[Pn -Pmf 



Pn ^ Pn + Pm 

II. n,'m 

X {{ll^n\dii^m){dj1prn\'4>n) + (V'n 1 9j V'm) (<9i-0m | V'n) ) (8) 

The inverse of the Pisher matrix provides a lower bound 
7 > H^^, on the covariance matrix 7^- = (XiXj) — 
{Xi){Xj) of global estimators of A, which is not generally 
achievable. On the other hand, the diagonal elements 
of the inverse Pisher matrix provide achievable bounds 
for the variances of single parameter estimators (at fixed 
value of the others) 



Var(A.) = 7.. > [H-'h 



(9) 



Let us now suppose to reparametrize the family of quan- 
tum states with a new set of parameters A = {Xj = 



Aj(A)}. We have = J^i^ijd^ with B^. 
and. in turn 



dXi /dXj 



Lj — ^ Bij Li 

i 

— E BirH{K)rsBj 



(10) 



i.e. H = BHB^ . 

In the following we address the problem of finding the 
bounds to the estimation of the entanglement between 
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the two subsystems A,B of a family of bipartite quan- 
tum states. The general strategy will be to start from 
the expression of the family of states in terms of a given 
set of " natural" parameters and then make a change of 
variable in order to write the state directly in terms of 
the chosen entanglement monotone A = e. Once this is 
achieved, the results obtained hold, at every fixed value 
of the entanglement, for the whole orbit of states that 
can be obtained by acting on the given one with local 
unitary operators Uab = Ua®Ub- Indeed, the latter, 
acting locally on the two subsystems A and B, they do 
not change the entanglement of the state; furthermore 
they do not change the value of the QFI. This follows 
from the observation that, if a given is the solution of 
([3]), the SLD that will correspond to — UAB0eU\g is 
given by = UabLeUab and, due to the cyclic prop- 
erty of the trace, F(e) = Tilg.L'^] = Tr[^eL?]. 

We finally notice that in order to assess the estimability 
of a given parameter the relevant figure of merit is given 
by the signal-to-noise ratio 

^(A) = < g(A) = X'H{X) , (11) 

rather than the variance itself. In particular the SNR 
of an estimator is relevant to assess its performances in 
estimating small values of the parameter. Eq. ([TT]) shows 
that the SNR i?(A) is bounded by the quantum signal- 
to-noise ratio (QSNR) Q{X) expressed in terms of the 
QFI. Upon taking into account repeated measurements 
we have that the number of measurements leading to a 
99.9% (Scr) confidence interval corresponds to a relative 
error 

2 _ 9Var(A) _ _9_ _ 9 

MA2 ^ MQ^ MX^H{\) 

Therefore, the number of measurements needed to 
achieve a 99.9% confidence interval with a relative error 
5 scales as 

^^(^) = Iq(a) 

In other words, a vanishing Q{\) implies a diverging num- 
ber of measurements to achieve a given relative error, 
whereas for a finite Q{X) the number of measurements is 
determined by the desired level of precision. In order to 
have a non- vanishing Q{\) for small value of the parame- 
ter A, the QFI should diverge at least as H{X) ^ A^^ for 
vanishing A. We notice that a similar quantity, namely 
XH{X), has been used in to asses estimation strategy for 
the parameter of a qubit depolarizing channel [s^ . 

III. TWO-QUBIT SYSTEMS 

In this section we analyze estimation of entanglement 
for families of two-qubit states and evaluate limits to pre- 
cision using the formalism developed in the previous sec- 
tion. At first we address the set of pure states and then 



consider families of entangled mixtures. In both cases we 
consider different measures of entanglement. 

A. Pure states 

We start by considering the set of pure states of two 
qubits. Upon exploiting the Schmidt decomposition 

^V^\0)a\0)b + VT^\1)a\1)b, (13) 

the whole family of pure states can be parametrized by 
a single parameter: the Schmidt coefficient q. Since for 
two-qubit pure states q is itself an entanglement mono- 
tone, all measures of entanglement can be expressed as 
a monotone function e ^ e{q). As a consequence, in 
order to determine the precision of estimation it suf- 
fices to evaluate the QFI H{q) and then use the rule for 
repametrization in Eq. pO|) . i.e. H{e) = H[q{e)][deq{e)]'^ . 
Since the states are pure, the SLD may be evaluated as 
Lq — 2dq\'i' q)('i/ q\; thc resulting Cramer- Rao bound and 
the QSNR read as follows 

Var(g) > H{q)-' = q{l - q) (14) 

Q{q) vanishes for vanishing q, thus indicating that any 
estimator of the Schmidt coefficient q becomes less and 
less precise for vanishing q. It is worth noting that the 
Schmidt coefficient q coincides with the only indepen- 
dent eigenvalue of the reduced density matrix qa(b) — 
diagjg, 1 — q\, which is diagonal in the Schmidt basis. 
Therefore the QSNR in Eq. ([15]) also imposes bound to 
the determination of the eigenvalue of qa{b)- Indeed, 
the same bound could have been obtained by applying 
the estimation machinery directly to Qa(b)- 

Let us now consider two different measures of entangle- 
ment for pure two-qubit states, i.e the negativity [33| 
and the (normalized) linear entropy eL = 2(1 — Tr[£i^]). 
In terms of the Schmidt coefficient q we have 

€^ = ,/r^ = 2^q{l-q) (16) 

We recall that the negativity is a good measure of entan- 
glement for generic two-qubit states, i.e it is an entan- 
glement monotone and it differs from zero iff the state 
is entangled, whereas the linear entropy is a good en- 
tanglement monotone only iff the state is pure. Upon 
expressing thc Schmidt coefficient as g = ^(1 — ^1 — e^) 
and using (fTO|) . we have Var(eN) > H{e^)~^ = 1 — e^, 
Var(eN) > H{ei^) = 4eL(l — Cl) and, in turn, 

Qie^)^-^'-^\l (17) 

Qi'-)- ^^'1,^) ^^-°e./4. (18) 

The optimal estimator for the Schmidt coefficient has a 
variance Var(g) which is minimum for g = 0, 1 (product 
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state) and maximum for (7=1/2 (Bell state) whereas 
for the two entanglement measure we have that Var(eN) 
is monotonically decreasing with e^; Var(eL) is minimum 
when the state is either in a product form (el — 0) or 
is maximally entangled (cl = 1) and is maximum in the 
"intermediate" case (el — 1/2). Despite the variances 
behave quite differently we have the same qualitative be- 
havior of the quantum signal-to-noise ratios and of the 
number of measurement necessary at fixed relative er- 
ror Mg. Indeed, in all cases the QSNR is an increas- 
ing function of the parameter and it diverges when the 
latter takes its maximum value (g = 1, = El = 1); 
Ms{e) diverges for e = = Cl = and then decreases 
monotonically, going to zero for the maximum value of 
entanglement e = = Cl = 1 ■ Moreover for vanishing 
entanglement the QNSR of the linear entropy estimator 
is vanishing slower than the corresponding quantity for 
the negativity. We conclude that the linear entropy is a 
more efficient entanglement estimator though, being the 
QSNR vanishing, the estimation is anyway inherently in- 
efficient. 

The above result can obviously generalized to the case 
of systems composed by a qubit and an TV-level system; 
indeed only two of the dimensions of the latter can be 
used to express the state: the reduced density matrices 
of both subsystems have only two non zero eigenvalues. 



new parameter {p,q) (/i,eN). The transfer matrix is 
given by 



B 



2^277^ 2^(2^.-1)2(2^-1-4) 




(22) 



2^(2m-1)(2m-1-4) 

and the inverse QFI matrix: 



The corresponding bound on the variance is thus given 
by: 



Var(eN) > (-H'(Ai, e^ 



1 



(24) 



which represents the bound to the precision of any entan- 
glement (negativity) estimation procedure performed at 
fixed purity ^. The result in Eq. p4|) is independent on 
the purity, no optimization procedure may be pursued, 
and it coincides with the bound obtained and discussed 
in the previous subsection for pure states. 

Let us now consider the family of Werner-like states 



Qpq = — ;— 1' 



(25) 



B. Entangled mixtures 

We now consider few families of mixed entangled states 
with different properties and show that they exhibit 
a common behavior concerning estimation of entangle- 
ment. The first family is described by the set of density 
matrices 



(19) 



where 



i.p=p|0,0)(0,0| + (l 
U{q) = exp{z arcos^/CTa; 



These states depend on two parameters {q^p) and are ob- 
tained via the action of the entangling operator U {q) on 
the classically correlated state Vp. Upon varying the pa- 
rameter p we may control the purity /i(p) = \ — 2p{l—p) 
of the state, while varying q we tune the amount of en- 
tanglement. The QFI matrix is diagonal with elements: 



H{p, q) = diag 



1 (1 - 2p)^ 



(20) 



The negativity of the state ([TO)) is given by the unique 
negative eigenvalue of the partially transposed state g^^ : 



e„ = 2Vg(l-<z)(l - 2p) 



(21) 



Upon inverting the above relation and the one for the 
purity, we reparametrize the set of states in terms of the 



obtained by depolarizing an entangled state \^q) of the 
form given in Eq. This set of states depends 

on two parameters p, q. As in the previous example 
upon varying the parameter p we may control the pu- 
rity ^{p) = (1 -I- 3p^)/4 of the state, while the amount of 
entanglement depends on both parameters. The eigen- 
values of Qpq depends only on p whereas the eigenvectors 
depends only on q. The QFI matrix is thus given by the 
diagonal form 



H{p, q) = diag 



l + (2-3p)p' q{l-q){l+p) 



(26) 



and the inverses of the diagonal elements correspond to 
the ultimate bounds to Var(p) and Yax[q) of any estima- 
tor of p and q, either at fixed value of the other param- 
eter or in a joint estimation procedure. Entanglement of 
Werner states may be evaluated in terms of negativity. 



max < 0, — 
^ 2 



1 r 



((1 + 4^9(9" 



1 



1 



(27) 



which implies that Werner states are entangled for (1 > 
p > [1 -f 4(7(1 — q)\~^ ■ Upon inverting Eq. (P7)) for p 
or q we may parametrize the Werner states using (e^, q) 
or (p, En) and evaluate the QFI matrices H{e^,q) and 
i/(p, Cn), their inverses and, in turn, the corresponding 
bounds to the precision of entanglement (negativity) esti- 
mation. The main results are that the ultimate bounds to 
the variance, and thus to the QSNR depend very slightly 
on the other free parameter (g or p). In other words, 
estimation procedures performed at fixed value of p or (7 
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respectively shows different precision, but the differences 
are negligible in the whole range of variations of the pa- 
rameters. We do not report here the analytic expres- 
sion of Q{e,^) at fixed p or q which is quite cumbersome. 
Rather, we show the behavior of Q(eN) in Fig. [1] On the 
left we show the QSNR Qg=o.5(eN) for q — 0.5, whereas 
on the right we show the ratio Qpie^) / Qq=o.5(^ti) for dif- 
ferent value of p (a similar behavior may be observed 
upon varying q). As it is apparent from the main panel 
the QSNR is a growing function of e^j vanishes for van- 
ishing negativity and diverges for maximally entangled 
states En = 1. The inset shows that there is almost no 
dependence on the actual value of p and q respectively 
and this prevents any possible optimization of the esti- 
mation procedure. For small we have Q{e^) ~ f{q)^^ 
and Q{e^) ~ g{p)e'^ respectively, where both the func- 
tions /(g) ~ 1 and ~ 1 are again very close to unit 
value for the whole ranges of variation of q and p. 



0.01 




FIG. 1: Quantum signal-to-noise ratio for the estimation of 
entanglement (negativity) of a two-qubit Werner state as a 
function of the negativity at fixed q — 0.5. The inset shows 
the ratio Qp(eN)/Q<j=o.5(eN) for different value of p. 

In all the cases we have considered, the QSNR is small 
for the most part of entanglement range and starts grow- 
ing only for highly entangled states. In other words, es- 
timation of entanglement is, on average, an inefhcient 
procedure. 



IV. TWO-QUTRIT BOUND ENTANGLED 
STATES 

In the previous section we have seen how the estima- 
tion of entanglement, as measured by negativity is a fairly 
inefficient procedure for weakly entangled states. Here 
we want to test how the QFI and the related bounds 
behave when one considers states that have inherently 
small amount of entanglement. A paradigmatic examples 
of such states are the so called bound entangled states, 
which exhibit non-classical correlations even if they sat- 
isfy the separability criterion based on partial transposi- 
tion of the density matrix [H, [13, HI] . The first example 



of bound entangled states is given by the following family 
of two spin-1 states 34 1 



Qa = 



l + 8a 



(Uo)(io| + ut)(iT l + |0i)(0i 



3a 



|OT>(OT| + |TO)(TO|) + — — 

1 



l + 8a 



|n>(n| 



(28) 



where 



\E)^^{\ ii) + |oo) + |TT)) 
|n> = A/^|Ti>^-/ 



tt). 



Since for all values of the parameter a Qa has a positive 
partial transpose (PPT), negativity cannot be used as 
a measure of the quantum correlations present in state. 
In order to estimate the entanglement we will use the 
scheme proposed in j39L [i^ . The latter is based on 
the following considerations. Given the sets of n non- 
commuting operators {Ai} and {Bi] acting locally on the 
subsystem A and B respectively one has a lower bound 
for the sum of the local uncertainties relations (LUR): 



^ 5Al > Ua and ^ SBf > Ub 



(29) 



where SOf ~ {Of) — {Oif' is the variance of the op- 
erator Oi. Since for all separable states on has that 
5{Ai + Bif' > Ua + Ub, the latter inequality set a nec- 
essary condition for a state to be entangled. The relative 
violation of the inequality defined as 



E.S{A, + B,Y 

Ua + Ub 



(30) 



can then be used as a measure of the quantum correla- 
tions present in the given state. The violation is neces- 
sary condition for the presence of the entanglement, thus, 
in order to effectively have and maximize such violation 
one can judiciously choose and optimize the choice of the 
sets {Ai} and {Bi}. The result of a possible optimization 
for the state ga is given in [40*1 and the corresponding rel- 
ative violation depends on the parameter a and is given 
by: 



3a^{l-a) 
4(2-ha)(l + 8a)2 



(31) 



The latter expression can be used to parameterize the 
state Qa in terms of ev{Qa) and then apply the QFI ma- 
chinery in order to obtain the desired bound on the esti- 
mation of relative violation of the LUR. The results are 
shown in Fig. [2l We first note that the relative violation 
of the LUR is small, i.e e^iQa) G [0,2/1125]. Nonethe- 
less, the number of measurements ^^(eu) is of the same 
order of those needed in the qubit case Ms{eK) in similar 
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conditions, and thus the overall efficiency of the estima- 
tion process is comparable for most part of the entangle- 
ment range. Finally, we notice that also for this family 
of qutrits the number of measurements Ms{eu) diverges 
for vanishing entanglement. 



MM 




0.0000 0.0005 0.0010 0.0015 



0.0000 0.0005 0.0010 0.0015 



FIG. 2; Number of measurement Ms(eu) needed to achieve 
a given relative error S = 10~^ (red), S — 10~^ (green), 
S = 10''^ (blue) in the the estimation of LUR entanglement 
measure of PPT state Qa in Eq.lgH}. Left: a £ [0,4/13], 
Right: a G [4/13, 1] 



V. GAUSSIAN STATES 

In this section analyze continuous variable systems and 
derive the bounds for the estimation of entanglement of 
two- mode Gaussian states 4l| . After a brief introduction 
we consider both pure states, i.e. twin-beam and the 
family of mixed states represented by squeezed thermal 
states (STS). Different measures of entanglement will be 
considered. 

The characteristic function of the state g is defined as 



is the 



x[g]iT) = Tr [^^(r)] where D{T) ^ exp 

displacement operator, defined in terms of the symplectic 
matrix 



k=l 



1 

-1 



(32) 



n 



and the vector R = (gi , pi , 
erators. — (ai, 6i, • • • , a 

coordinates and we have [Rk,Rh\ = ^^kh- Two- mode 
Gaussian state are those with a characteristic function of 
the form 



7 Qn 1 Pn 

) of canonical op- 
■n,bn) denotes the Cartesian 



x[e](r) = exp(-ir^ar 



(33) 



into the standard block form : a 



where A ■ 



A C 
C B 

Diag(a, a), B = Diag(6, 5), and C — Diag(c+,c_). For 
two-mode Gaussian states PPT condition is necessary 
and sufficient for separability and thus the entanglement 
properties of the state are encoded in the least symplectic 
eigenvalue [i^j 



d- = \/ {a - c+)(a + c_) 



(35) 



where the symplectic spectrum of a^^ can be evaluated 
by finding the eigenvalues of ^~^a^^ . In this framework 
the PPT criterion can be cast in term of the smallest sym- 
plectic eigenvalue i.e., p is separable iff d_ > 1/2. Indeed, 
(i_ is itself an entanglement monotone. Furthermore, all 
the different entanglement measures for symmetric Gaus- 
sian states that have been proposed [H, [H, turns 
out to be a monotone function of the smallest symplectic 
eigenvalue. Since we are interested in the estimation of 
entanglement, we may first study the estimation of the 
symplectic eigenvalue d- and then use repametrization 
in order to asses the performances of other entanglement 
monotones. In particular, we will focus on the following 
measures: 



en[d^) = max{0, — \\i2d-} 
eL(d-, 
es(d-] 
e^{d-] 



1 

l + Adl. 
1 - 2d 

(1 - ^/2df 
l + 2d 



(36) 
(37) 
(38) 
(39) 



where 



where the expressions for the linear entropy El has been 
obtained in the pure state case. In particular El ^md the 
logarithmic negativity will be used for pure states, 
whereas the Bures distance-based measure and eg [i^ 
will be used for mixed states. Notice that Eb is a good 
measure of entanglement just for symmetric two-mode 
Gaussian states and that Eq. ((38)) has been obtained for 
this particular class of states. In order to evaluate the 
QFI, one has first to determine the actual expression of 
d_. Then, one expresses the elements of the covariance 
matrix in terms of the chosen entanglement monotone 
and proceeds with the repametrization rules described in 
Section nil 



Pure states 



^kh = -{{Rk,Rh}) - {Rk){Rh) 



(34) 



is the covariance matrix of the second moments and X 
denotes mean values. The second term in the exponential 
does not contain any information about entanglement, 
and can be set to zero via local operation. The covariance 
matrix completely characterize the state and by means 
of local symplectic transformations can be transformed 



Here we address estimation of the entanglement of pure 
two modes Gaussian states, i.e. twin-beam. These are 
defined by the following relations between the elements 
of the covariance matrix: a = b i.e., the states are sym- 
metric and c+ = — c_ = -y/a^ — 1/4. The states can thus 
be described by the single parameter a or, by inverting 
d_ ~ a — ^ a} — 1/4 they can be completely described 
by their entanglement content. For Gaussian states the 
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evaluation of the SLD and the QFI, besides the use of 
Eqs. (O, may be pursued using phase-space techniques. 
In fact, for pure states we have = 2d^Q^ and this al- 
lows to directly evaluate the characteristic function of the 
SLD as follows 



(40) 



where = [T^^ , Ts^). The corresponding QFI H{e) = 
Ty[qL1] is given by 



j(e,ri,r2) (41) 



27r 27r 

where the integrand function reads as follows 

I = x[Wi)x[L,]{T2) Tr [bD\T^)D\T2)] (42) 

We now use the relations -D[r]^ = T] and 

D{Ti)D{T2) = D{Ti+T2)g{Ti,T2) and we rewrite gS]) 

as 

1 = {r'^d.a.Ti) {r^d,a,T2) ■ 

x,.(ri)xe,(r2)xe,:(-ri-r2)g(ri,r2) 
= (r^Sir) (r^Sar) exp(^-ir^Ar^ (43) 



where we have introduced the matrices; 



S = 



2(7 a 
a 2(7 



.Si = 



d,a 




,S2 = 





acr 



and where A = S -I- iT, with T = ^ CTj, ® I (g) cr^ in 
terms of Pauli matrices. The result of the integration QFI 
is a function of £,d^a^ and can be expressed in various 
ways depending on the entanglement monotone that one 
chooses to estimate. By setting a = ^ cosh Cn in the 
covariance matrix, one can use the logarithmic negativity 

. In this case one finds that the QFI is independent on 
the entanglement content of the state. 

If one uses directly the symplectic eigenvalue the 
QFI now depends on the entanglement monotone: 



H{d^) = dZ^ 



(44) 



and the minimal variance in the estimation can be ob- 
tained in the limit of infinite entanglement i.e., d- — 0. 
Moreover we observe that, while for the logarithmic neg- 
ativity one sees that the QSNR Q{etq) is simply propor- 
tional to e^, if we consider the least symplectic eigenvalue 
we have Q{d-) = 1 over the whole range of variation, i.e. 
the estimation procedure can be done efficiently either for 
highly entangled states and weakly entangled ones. 

As a matter of fact, twin-beam may be also written in 
the Fock basis as l^*) = /n|'^)|'^) where, in terms of 
the log-negativity or the linear entropy one may write 



J (cosh En — 1)" 

'(cosheN + 1)1+") 



\l+n 



(2 



(45) 



Using this representation we may directly exploit Eqs. 
([7]): for a generic parameter x and fn{x) G M wc have 
(^fja^r*) = and thus 



Hix) = - ^(5,/„(x))2 



(46) 



Using the above equation one recover the result for the 
log-negativity and may evaluate the QFI in terms of the 
linear entropy obtaining 



H{e^) = (4(2-eL)(eL-l)'eL 



B. Entangled mixtures 



(47) 



We now analyze entanglement estimation for a rele- 
vant family of mixed Gaussian states labeled by two inde- 
pendent parameters. The symmetric two-mode squeezed 
thermal states (STS) are given by 



0ST 



Suir, ^){iyN, ® mJSl2{r, cj,) (48) 



and represents a two parameter family qst — QsT{Nt,r) 
obtained from symmetric two-mode thermal state with 
Nt thermal photons for each mode by the action of the 
two-mode squeezing operator Si2{r) = exp{r(a|a2 — 
0102)}. The family in Eq. occurs when one consid- 
ers the propagation of twin-beam in a noisy channel or 
the generation of entanglement from a noisy background 
[HI, and represents the CV generalization of the fam- 
ily of entangled mixed states introduced in Eq. p9p . 
We evaluate the corresponding Fisher information ma- 
trix H{r, Nt) and obtain a diagonal matrix 



H{r,Nt) = diag{8- 



l + 2Nt{l + Nt)' Nt{l + Nt)J ' 

(49) 

Let us now consider the smallest symplectic eigenvalue 
d-{r, Nt) and the purity of the state fi{Nt) 



-{l + 2Nt) M 



1 



2Nt + l 



Upon inverting the above equations we may 
reparametrize the set of states in terms of the new 
parameters (d_,/i), the transfer matrix B being given 

by 



B 



The new QFI matrix H{d^,^) is calculated by means 
of Eq. (jlOp and the bound on the covariance matrix 
7(d_,/i) > H{d^, fj,)^^ is established by its inverse 




Hid^,^i)-' 



J2 c'-A'(l-M^) 

^(i-zi^) 



(50) 
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The lower bounds on the variance for the symplectic 
eigenvalue is given by 



(51) 



and represents the limit to the precision of any estimator 
of d- at fixed purity fi. In particular, we observe that this 
bound does not depend on the purity, and coincides with 
the bound in Eq. (|44|) obtained for pure states. Therefore 
also for this class of states the QSNR is Q{d-) = 1 and 
hence d- can be always estimated efficiently. 

Let us now consider a generic measure of entanglement 
e = e((i_). Upon using Eq. ([TU)l we may show that the 
reparametrization — s- (e,^) leads to 



(52) 



dd 



Var(e) > (JJ(e, = (^^ j (Jf(d_, m)-')ii 

(53) 

Let us consider the two monotone functions of the sym- 
plectic eigenvalue es(d-) and e^id-) introduced in Eqs. 
([55)1 and ([55]) . The symplectic eigenvalue can be ex- 
pressed in terms of the measures as 



dies) 
d(e^) 

thus leading to 



1 -e« 



2(l-eB)2 



Var(es) > (1 - e^ f 
Var(e.) > ^^i' ' ' 

We notice that Var(es) and Var(eB) show different be- 
havior; in particular, while the bound on Var(es) van- 
ishes only when es is maximum (es — 1), the bound on 
Var(eB) reaches zero both when Eb is maximum (eB = 1) 
and when is minimum (eB = 0) and presents a maximum 
for eB = 1 - 1/V2. 

We finally evaluate the QSNR for the measures of en- 
tanglement introduced, obtaining 



Q{es) < 



Q{e 



< 



4ep 



(54) 



«B^0 



(l-eB)2(2-eB) 



2ep 



The two QSNRs are increasing function of entanglement, 
vanish for zero entanglement and diverge for maximally 
entangled states. In turn, the numbers of measurements 
Ms{e b) and Ms{es) vanish for maximum entanglement 
and diverge for vanishing entanglement. The QSNR of 
Cb is vanishing slower than the corresponding quantity 



for eg and therefore we conclude that the measure based 
on the Bures distance is more efficiently estimable com- 
pared to the linear measure es . On the other hand, being 
the QSNR vanishing, the estimation is anyway inherently 
inefficient. 



VI. CONCLUSIONS 

Entanglement of quantum states is not an observable 
quantity. On the other hand, the amount of entangle- 
ment can be indirectly inferred by an estimation proce- 
dure, i.e. by measuring some proper observable and then 
processing the outcomes by a suitable estimators. In this 
paper we have established a first approach to the estima- 
tion of the entanglement content of a quantum state and 
to the search of optimal quantum estimators, i.e those 
with minimum variance. Our approach is based on the 
theory of local quantum estimation and allows, upon the 
evaluation of the quantum Fisher information, to derive 
the ultimate bounds to precision imposed by quantum 
mechanics. We have applied our analysis to several fam- 
ilies of quantum states either describing finite size sys- 
tems or continuous variables ones, and have considered 
different measures in order to quantify the amount of en- 
tanglement. 

For the case two-qubit pure state we have found that 
any procedure to estimate entanglement (either quan- 
tified by negativity or by linear entropy) is efficient 
only for maximally or near maximally entangled states, 
whereas it becomes inherently inefficient for weakly en- 
tangled states. In particular, the number of measure- 
ments needed to achieve a 99.9% confidence interval with- 
ing a given relative error diverges as far as the value of 
entanglement becomes small. The same results hold also 
for families of mixed states, remarkably for the orbit of an 
entangling unitary an for a general class of Werner-like 
states. Indeed in all the examples we have considered the 
presence of other free parameters besides entanglement, 
though changing the QFI, does not affect the estimation 
precision, i.e. the value of the relevant element of the 
inverse QFI matrix. In turn, this also prevents the pos- 
sibility of further optimizing the estimation procedure. 

On the other hand, we have showed that for an impor- 
tant class of states whose entanglement of distillation is 
zero (PPT bound entangled states), the use of an opti- 
mized measure of quantum correlation i.e., the relative 
violation of local unitary relations introduced in , re- 
sults in a more efficient estimation procedure, with preci- 
sion comparable with those achievable in the estimation 
of entanglement through negativity. 

In the case of continuous variable Gaussian states we 
have shown that the estimation of the least symplec- 
tic eigenvalue d_ of the covariance matrix may be per- 
formed with arbitrary precision at fixed number of mea- 
surements, independently on the value of d_ itself and for 
both pure states and mixed states. If we rather introduce 
other measures of entanglement proposed in literature, in 
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particular the logarithmic negativity for pure states and 
the one based on the Bures distance |43| for the symmet- 
ric squeezed thermal (mixed) states, we observe the same 
behavior obtained in the discrete variable case: the esti- 
mation is efficient only for maximally entangled state and 
inherently inefficient for weakly entangled states. There- 
fore it is apparent that for continuous variable systems, 
the efficiency of the estimation strongly depends on the 
measure one decides to adopt. 

In conclusion, upon exploiting the geometric theory 
of quantum estimation we have quantitatively evaluated 
the ultimate bounds posed by quantum mechanics to the 
precision of entanglement estimation for several families 
of quantum states. To this aim we used the quantum 
Cramer-Rao theorem and the explicit evaluation of the 
quantum Fisher information matrix. We have also given 
a recipe to build the observable achieving the ultimate 
precision in terms of the symmetric logarithmic deriva- 
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